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^\ ' Abstract. Observations of the radio galaxy MRC 0316-257 at z = 3.13 and the surrounding field are presented. Using narrow - 

, and broad-band imaging obtained with the VLT, 77 candidate Lya emitters with a rest-frame equivalent width of > 15 A were 

selected in a ~ 7'x7' field around the radio galaxy. Spectroscopy of 40 candidate emitters resulted in the discovery of 33 
emission line galaxies of which 31 are Lya emitters with redshifts similar to that of the radio galaxy, while the remaining 
two galaxies turned out to be [On] emitters. The Lya profiles have widths (FWHM) in the range of 120-800 kms~', with 
a median of 260 kms~'. Where the signal-to-noise was large enough, the Lyor profiles were found to be asymmetric, with 
apparent absorption troughs blueward of the profile peaks, indicative of absorption along the line of sight of an H i mass of at 
least 2 x 10 2 - 5 x 10 4 Mq. Besides that of the radio galaxy and one of the emitters that is a QSO, the continuum of the emitters 
is faint, with luminosities ranging from 1.3 L* to < 0.03 L* . The colors of the confirmed emitters are, on average, very blue. 

The median UV continuum slope is/3 = -1.76, bluer than the average slope of LBGs with Lya emission (/? 1.09). A large 

fraction of the confirmed emitters (~ 2/3) have colors consistent with that of dust-free starburst galaxies. Observations with 
the Advanced Camera for Surveys on the Hubble Space Telescope show that the emitters that were detected in the ACS image 
have a range of different morphologies. Four Lya emitters (~ 25%) were unresolved with upper limits on their half light radii 
of r h < 0.6 - 1.3 kpc, three objects (~ 19%) show multiple clumps of emission, as does the radio galaxy, and the rest (~ 56%) 
are single, resolved objects with r h < 1.5 kpc. A comparison with the sizes of Lyman break galaxies at z ~ 3 suggests that the 
Lya emitters are on average smaller than LBGs. The average star formation rate of the Lya emitters is 2.6 Mq yr~' as measured 
5^ ' by the Lya emission line or < 3.9 Mq yr~' as measured by the UV continuum. The properties of the Lya galaxies (faint, blue 

and small) are consistent with young star forming galaxies which are still nearly dust free. 

The volume density of Lya emitting galaxies in the field around MRC 0316-257 is a factor of 3.3^ larger compared with 
the density of field Lya emitters at that redshift. The velocity distribution of the spectroscopically confirmed emitters has a 
dispersion of 640 kms~', corresponding to a FWHM of 1510 kms~', which is substantially smaller than the width of the 
narrow-band filter (FWHM ~ 3500 km s ). The peak of the velocity distribution is located within 200 km s of the redshift of 
the radio galaxy. We conclude that the confirmed Lya emitters are members of a protocluster of galaxies at z ~ 3.13. The size 
of the protocluster is unconstrained and is larger than 3.3 x 3.3 Mpc 2 . The mass of this structure is estimated to be > 3-6 xlO 14 
Mq and could be the progenitor of a cluster of galaxies similar to e.g. the Virgo cluster. 

Key words. Galaxies: active — Galaxies: high-redshift — Galaxies: evolution — Galaxies: clusters: general — Cosmology: 
observations — Cosmology: early Universe 
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* Based on observations carried out at the European Southern 
Observatory, Paranal, Chile, programs LP167.A-0409 and 68.B-0295. 



1. Introduction 

Within Cold Dark Matter (CDM) scenarios the first stars and 
stellar systems form through gravitational infall of primordial 
gas in large CDM halos (e.g. lWhite & Reeslfl978i) . Numerical 
simulations suggest that as these halos merge they form vast, 
web-like networks of young galaxies and ionized gas (e.g. 
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Bau gh etail Il998l) . The most massive galaxies, and the rich- 
est clusters emerge from regions with the largest overdensities. 
Although clusters of galaxies have been studied extensively out 
to z ~ 1.3 (e gjRosatietallll999tlpehaCecaeta"lll200Q| 
Stanford et all 120021 iBlakeslee et all l2003bt iMauehan et all 
20031 iToftet all 12004 . the epoch of cluster formation is still 



an open question due to the difficulty in identifying their pro- 
genitors in the early Universe. 

During the last decade, evidence has mounted that the 
most powerful high redshift radio galaxies (HzRGs; z > 
2) are progenitors of brightest cluster galaxies and are lo- 
cated in dense environments. HzRGs are amongsTthe bright- 
est and presumably most massive galaxies dJarvisetallboOll 
iDe Breuck e7alll2002tlzirm et allkooH) . They have high star 
formation rates (> 100 M» yr' 1 ), bas ed on deep spectra of their 
UV continuum ("e.g.lDev et allll997l) and the detections of dust 



(e.g.lArchibald et alll200ltlStevens et a"Hl2003HReuland et all 



2004 |) and extended CO e mission J^padopoul osetaiT BoOO ; 
De Breuck et all l2003a b). Furthermore, radio galaxies at red- 
shifts between 0.5 and 1.5 are known to predom i nantly 
lie in moderat ely rich clusters (Hill & Lilly, Il99lt iBestL 
l2000t iBest etail 120031) . At higher redshifts (z > 2), some 
radio galaxies were found to posses companion galaxies 
jLe Fevre et allll996tlPascarelle et allll996tlRottgering et all 
ll996HKeeletall ll999l) . Also. 20% of the HzRGs have ex- 
treme radio rotation measures (> 1000 radrrT 2 ), giving an in- 
dication that these radio galaxies are surrounded by dense hot 
gas JCarilli et all 119971 lAthreva et all Il998l IPentericci et all 
l2000bl) . 

To search for direct evidence of the association of a clus- 
ter or a forming cluster (protocluster) with a radio galaxy, we 
conducted a pilot project on the Very Large Telescope (VLT) 
aimed at finding an excess of Lya emitters around the clumpy 
radio galaxy PKS 1 138-262 at z = 2.16. Narrow-band imaging 
resulted in a list of ~ 50 candidate Lya emitters (IKurk et all 
l2000U2004 . Subsequent multi-object spectroscopy confirmed 
14 Lya emitting galaxies and one QSO whos e velocities were 
within 1000 km s" 1 of the central radio galaxy JPentericci et all 
2000al lKurketalll2004 . The volume density of Lya emitters 
near PKS 1138-262 was found to b e a factor 4 . 4 + 1 .2 times 
that of Lya emitters in blank fields dKurk et ail 12004). Using 
near-infrared narrow- and broad-band images of the field, sig- 
nificant populations of Ha emitters at the redshift of the radio 
galaxy and extremely red objects were found. Also, Chandra 
observations revealed a n excess of soft X-ray sources in the 
field of PKS 1 138-262 JPentericci et alll2002h . indicating that 
several AGN are present in the protocluster. 

As shown by the study of the overdense region near PKS 
1138-262, distant protoclusters provide ideal laboratories for 
tracing the development of large scale structure and galaxy evo- 
lution. To further study the formation of large scale structure in 
the early Universe and to investigate the evolution of galax- 
ies in dense environments, we initiated a large program on the 
VLT to search for Lya emitting galax ies around luminous ra - 
dio galaxies with redshifts 2 < z < 5 ( Ven emans et all 12003). 
The goals were to find protoclusters of galaxies, determine the 
fraction of HzRGs associated with protoclusters and study the 
properties of protoclusters and their galaxies. The first result 



was the discovery of a protoclust er around the r adio galaxy TN 
J1338-1942 at z = 4.1 ( I Venemans et alll2002h . Deep imaging 
and spectroscopy revealed 20 Lya emitters within a projected 
distance of 1.3 Mpc and 600 kms -1 of the radio galaxy. By 
comparing the density of Lya- emitters in the protocluster to 
the field, the galaxy overdensity was claimed to be 4.0 ±1.4 
and the mass of the s tructure was estimated to be ~ 10 15 Mo 
JVenemans et alll2002l) . 







Here we report on observations of the radio galaxy MRC 
0316-257. This 1.5 Jy radio sour ce was listed in the 408 MHz 
Molonglo Referen ce Catalogue jLarge et all Il98lh and opti- 
cally identified by McCarth y et all Jl990i>. Its discovery spec- 
trum yielded a redshift of 3.13 [McCarth y et all Il990) . This 
object was included in our program because the redshift of 
the radio galaxy shifted the Lya line into one of the narrow- 
band imaging filters available at the VLT. Also, it already 
had two spectroscopically confirmed Lya emitting compan- 
ions (Le Fevre et all 1 199 6. hereafter LF96), an indication that 
the radio galaxy is located in a dense environment. Further, 
the redshift of the radio galaxy of 3.13 allows for an efficient 
search for Lyman Break Galaxies (LBGs) and for [O ra] /15007 
A emitters using a K-bwd narrow-band filter, which is avail- 
able in the Infrare d Spectrometer and Array Camera (ISAAC, 
lMoorwoodlll997l) at the VLT. 

Besides observing MRC 0316-257 with the VLT as part 
of our large program, we made additional observati ons of the 
field w ith the Advanced Camera for Surveys (ACS: lFord et all 
Il998l) on the Hubble Space Telescope (HST) 1 to study the sizes 
and morphologies of the detected galaxies. 

This paper is organized in the following way. In Sect. |2 
the imaging observations and data reduction are described and 
Sect.[3]discusses how candidate Lya emitters in the field are de- 
tected. The spectroscopic observations and the results are pre- 
sented in Sect. |4] The properties of the Lya emitters are ana- 
lyzed in Sect. [5] and details of individual emitters are presented 
in Sect. [6] Evidence for the presence of a protocluster in the 
field is discussed in Sect.0 and the properties are presented in 
Sect. [8] In Sect. [9]the nature of the Lya emitters is discussed, 
followed by a description of the implications of a protocluster 
atz= 3.13 in Sect.Ql 

T hrough out this article, magnitudes are in the AB system 

Jokelll974 . using the transformations Vab — VVega + 0.01 and 
Iab = ^Vega + 0.39 (Bessell, 1979). A A-dominated cosmology 
with H = 65 kms" 1 Mpc~', Q M = 0.3, and Q A = 0.7 is 
assumed. In this cosmology, the luminosity distance of MRC 
0316-257 is 28.8 Gpc and 1" corresponds to 8.19 kpc at z = 
3.13. 



1 Based on observations made with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, which 
is operated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS 5-26555. These obser- 
vations are associated with program #8183 
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Date 


Telescope 


Instrument 


Mode 


Optical element 


Seeing 


Exposure time 


2001, September 20 and 21 


VLT UT4 


FORS2 


Imasins 


Bessel V 


077 


4 860 s 


2001, September 20 and 21 


VLT UT4 


FORS2 


Imaging 


OIII/3000 


077 


23 400 s 


2001, September 22 


VLT UT4 


FORS2 


MOS fl 


GRISJ400V 


~ 1'.'5 


12 600 s 


2001, October 18 


VLT UT4 


FORS2 


MXU*, mask I 


GRISJ400V 


l'.'O 


10 800 s 


2001, October 18, 19 and 20 


VLT UT4 


FORS2 


MXU*, mask II 


GRISJ400V 


l'.'O 


29 100 s 


2001, November 15 and 16 


VLT UT3 


FORS1 


PMOS c 


GRIS.300V 


078 


19 800 s 


2002, July 18 


HST 


ACS 


Imaging 


F814W 




6300 s 


2002, September 6, 7 and 8 


VLT UT4 


FORS2 


Imaging 


Bessel / 


a' 7 


4 680 s 



" Multi-object spectroscopy mode, performed with 19 movable slitlets with lengths of 20"-22". 

Multi-object spectroscopy mode with a user-prepared mask. 
c Spectropolarimetry mode using 9 movable slitlets of 20". 



2. Imaging observations and data reduction 

2.1. VLT imaging 

An overview of the observations is shown in Table 1. On 2001, 
September 20 and 21, narrow- and broad-band imaging was 
carried out with the 8.2 m Yepun (VLT UT4) to search for 
Lya emitting galaxies around MRC 0316-257. The instrument 
used was the FOcal Reducer/low dispers ion Spectrograph 2 
(FORS2; App enzeller & Rupprechtl [l992) in imaging mode. 
For the narrow -band imaging the OIII/3000 filter was used with 
a central wavelength of 5045 A and full width half maximum 
(FWHM) of 59 A, which samples the Lyq line from the radio 
galaxy, which is redshifted to 5021 A (McCarthy et al., 1990, 
LF96). To measure the UV continuum near the Lyo- line, the 
field was imaged with broad-band filter Bessel V with a cen- 
tral wavelength of 5540 A and FWHM of 1 1 15 A. The detector 
was a SiTE CCD with 2048x2048 pixels. The pixel scale was 
072 per pixel, resulting in a field of view of 6'8x6'8. A year 
later, on 2002, September 6, 7 and 8, broad-band images of the 
field were taken in the Bessel / filter, with a central wavelength 
of 7680 A and a FWHM of 1380 A. The instrument was again 
FORS2, but the detector was replaced by two MIT CCDs with 
2048x2048 pixels each. The gap between the two CCDs was 
~ 4 arcsec. The pixel scale of the MIT CCDs was 07125 per 
pixel. To decrease the readout time, the pixels were binned by 
2x2, resulting in a spatial scale of 0725 pixeL 1 . The field 
of view was restricted by the geometry of the Multi-Object 
Spectroscopy (MOS) unit, and was 6'8x6'8. 

The observations in the narrow-band were split into 13 sep- 
arate exposures of 1800 s, in V-band into 27 separate 180 s 
exposures and in the /-band 26 exposures of 180 s were taken. 
The individual exposures were shifted by ~ 15" with respect 
to each other to facilitate identifying cosmic rays and remov- 
ing residual fiat-field errors. 

All nights except for 2001, September 20 were photomet- 
ric, and the average seeing was 0765 - 077 in the narrow-band, 
V and / images (see Table 1). For the flux calib ration, the 
spectrophotometric standard sta r LTT 1788 JStone & Baldwin, 
1 19831 iBaldwin & StoneL 1 1984 was observed in th e V-band 
and th e photometric standard stars in the field SA98 JLandold 
1 19921) were used to calibrate the /-band images. 



2.2. Data reduction of VLT data 

The VLT images were reduced using standard routines within 
the reduction software package IRAF 2 . The reduction steps in- 
cluded bias subtraction, flat fielding using twilight sky flats and 
illumination correction using the unregistered science frames. 

The magnitude zero-points derived from different standard 
stars were consistent with each other within 0.02 magnitude. 
To derive the zero-point of the narrow-band image, the magni- 
tude of the ~ 400 brightest objects in the field were measured 
in the V and /-band images. These objects had a signal-to-noise 
of at least 25 in both V and /-band images. A magnitude limit 
of mi > 20 was set to reject saturated stars. Narrow-band mag- 
nitudes were derived from the V and /-band magnitudes as- 
suming a powerlaw spectral energy distribution for these 400 
objects. With these derived narrow-band magnitudes and the 
associated counts in the narrow-band image, the zero-point of 
the narrow-band image was computed. The rms of the com- 
puted zero-point was 0.006 mag. 

All magnitudes we re corrected for gal actic extinction 
which was estimated by Schlege l'et alJ Jl99 8) to have a value 
of E(B — V) = 0.014 mag. The measured 1 <x limiting magni- 
tudes per square arcsecond were 28.35, 28.90 and 28.69 for the 
narrow-band, V-band and /-band respectively. 

Astrometric calibr ation was performed us ing the USNO- 
A2.0 catalog jMonet et allll998llMonelll998l) from which 20 
stars were identified in the field. This resulted in a fit with a 
typical error in the right ascension and declination of 07 17. The 
astrometric accuracy of the images is dom inated by the u ncer- 
tainty of the USNO-A2.0 catalog of 0725 jDeutschLll999i) . The 
VLT images were registered in the following way. Because the 
V-band and narrow-band images were taken with the same de- 
tector, a simple pixel shift was sufficient to align the images, 
using the positions of a few stars over the field. The /-band 
frames were taken with the MIT CCDs, which had a differ- 
ent pixel scale (0725 pixeL 1 ) as compared to the SiTE CCD 
(072 pixeL 1 ). The distortions were also different, and rescal- 
ing the /-band image to the same pixel scale as the V-band 



2 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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image resulted in positional errors up to a few arcseconds in 
the corners of the /-band image. Instead, the position of a few 
hundred objects detected in both the /-band and V-band images 
with a signal-to-noise greater than 15, were used for the align- 
ment. This way the positional error of objects in the /-band im- 
age dropped to 0704. S ubsequently, the algorithm DRIZZLE 
jFruchter & HookLl2002i) was used to map the /, V and narrow- 
band images on new frames with a common pixel scale of 07 16. 

The area of the reduced images was 46.7 arcmin 2 . Due 
to the presence of two bright stars in the field, the area that 
could be used for detecting candidate Lya emitters was 45.75 
arcmin 2 . The width of the narrow-band filter in redshift is 0.049 
at z ~ 3.13 and the volume probed by the filter at z = 3.13 is 
9331 Mpc 3 . 
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2.3. Hubble Space Telescope imaging and reduction 

A part of the field imaged by the VLT was observed in 2002 
July with the ACS on board the HST as part of an imaging 
program of HzRGs. The 3.'4x3?4 field of view of the ACS was 
chosen to include not only the radio galaxy but also as many 
confirmed Lya emitters as possible (see Fig.ll4lfor the position 
of the ACS field within the FORS field). The field was imaged 
in the F814W filter (hereafter /814) with a central wavelength 
of 8333 A and a width of 251 1 A. The total exposure time was 
630 s. The images wer e reduced using the ACS GTO pipeline 
( Blakeslee et alll2003ah. 



3. Detection and selection of candidate emitters 

3.1. Source detection 

For the detection and photometry of objects in the images, the 
program SExtractor (version 2.2.2, Bertin & Arnoutsl [1996) 
was applied. The narrow -band image was used to detect the ob- 
jects. Because some of the Lya emitters remained undetected in 
the broad-band images (see Table 4), this was preferred above 
a combination of the narrow-band and broad-band images as 
dete ction image, w hich is favoured by some other groups (e.g. 
iFvnbo et all Bo02). Detected objects in the narrow-band image 
were defined to have at least 15 connected pixels with values 
larger than the rms sky noise. This resulted in a list of 3505 ob- 
jects, of which 3209 had a signal-to-noise greater than five. To 
assess the completeness of the catalog, artificial and real galax- 
ies were added to the narrow-band image and recovered. The 
galaxies had various sizes, the smallest galaxies had a half light 
radius r/, ~ 0.4", similar to that of stars in the field (unresolved 
objects), the largest galaxies had a half light radius r/, ~ 0.9", 
roughly 2.5 times that of stars. We found that the completeness 
depended on the size of the galaxies that were added to the 
image as shown in Fig. ^ The limit where half of the galax- 
ies were recovered ranged from a magnitude of m n b < 26.25 
for unresolved objects to m n \, < 25.25 for the largest objects. 
The 90% recovery limit was m n b ~ 26.0 for unresolved and 
m n b ~ 25.0 for the largest objects. 



Fig. 1. Fraction of galaxies recovered in the narrow-band image 
versus magnitude for various galaxy sizes. 

3.2. Photometry 

To measure the flux of the detected objects, the double image 
mode of SExtractor was employed. In this mode, SExtractor 
detects objects in one image, and carries out the photometry on 
a second image. In our case, the narrow-band image was used 
for the detection of the objects and the photometry was done 
on the narrow-band and broad-band images. Of each detected 
object in the narrow-band the flux was measured in two aper- 
tures: a circular aperture to compute the colors of the object, 
and an elliptical aperture to estimate the total brightness of the 
object. The radius of the circular aperture (/? apel ) depended on 
the isophotal area of the object (A; so ), which is the area of pixels 
with values above the rms sky noise and is an output parameter 
of SExtractor: /? aper = y]A lso ln. A minimum aperture radius of 
07525, 1.5 times the radius of the seeing disc, was set to avoid 
very small apertures. The maximum radius was set to 4 times 
the radius of the seeing disc to avoid overlapping apertures due 
to neighbouring galaxies. The shape and size of the elliptical 
aperture was derived from the object's light distribution. The 
ellipticity s and position angle of the object were computed 
from the second order moment of the light distribution. Using 
the first moment r\ , the elliptical apert ure had ma jor and minor 
axes of kri/s and ekr\ (Ber tin & Arno uts. 1996). The scaling 
factor k determines the size of aperture and is a free parame- 
ter in SExtractor. We tested SExtractor on a set of images and 
compared SExtractor's output magnitudes for values of k in the 
range 1.0 < k < 2.75. It was found that a scaling factor of 
k =s 1 .75 both optimized the signal-to-noise and minimized the 
fraction of the flux of the object outside the aperture. 

The aperture used to measure the total flux of an object was 
the elliptical aperture, except when more than 10% of the pixels 
in the elliptical aperture was significantly effected by bright and 
close neighbours (SExtractor output parameter FLAGS equals 
1) or when the object was originally blended with another one 
(FLAGS equals 2). In those cases the circular aperture was used 
to derive the total flux. 

To estimate the fraction of the total flux of the object falling 
outside the (elliptical) aperture, Monte Carlo simulations were 
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Fig. 2. Difference between measured magnitude and input mag- 
nitude as a function of input magnitude for objects with various 
sizes. The difference is larger for fainter and/or larger objects. 

performed. Galaxies in the magnitude range 21 < m n b < 28 
with various shapes (Gaussian and elliptical profiles) and sizes 
(half light radii 0.4" < r;, < 0.9") were added to the narrow- 
band image and recovered. The number of galaxies added to 
the image was limited to 40 to avoid overcrowding. This rou- 
tine was repeated until 8000 galaxies per magnitude bin were 
simulated. By comparing the measured magnitudes to the in- 
put magnitude of the simulated galaxies, the fraction of the flux 
that was outside the aperture could be estimated to correct the 
aperture magnitudes. The simulations showed that this fraction 
depends on both the original size of object and the magnitude 
of the object: the fraction of the flux outside the aperture is 
higher for a large and/or faint object, compared to a compact 
and/or bright object (see Fig.[5J. At faint magnitudes the cor- 
rection becomes smaller again, because these objects have a 
larger probability to be detected if they coincide with a peak in 
the noise. For the brightest objects the fraction of the flux out- 
side the aperture is constant at a value of ~ 11%. To avoid an 
overestimation of the magnitude correction, we decided to use 
only point sources to measure the correction (see Fig.|2j. The 
magnitude correction applied to the sources in our field was 
< 0.1 for objects with m n b < 21, rising to ~ 0.25 for objects 
with m n b ~ 26. 

3.3. Selection of candidate Lya emitters 

An efficient method of detecting Lya emitting galaxies at high 
reds hift is to selec t objec ts with a large line equivalent width 
(e.g. Co wie & Hul Il998l) using the narrow- and broad-band 
photometry. The observed equivalent width £W b s of a Lya line 
is defined as: 



(1) 



with /*Lya the flux of the Lya line and C\^ a the UV continuum 
at the wavelength of the Lya line. Assume that a Lya line is 
observed in both a narrow-band filter and a broad-band filter, 
then the flux density in the narrow-band (/i, n b) and broad-band 
(/i,bb) can be described as: 



Anb - C(A - /ieff.nb) + F Lya /AA n \, 
Abb = C{A = /i e ff,bb) + ^LyoV A/lbb, 



(2) 
(3) 



where C is the UV continuum and A/lbb(nb) is the width of 
broad-band (narrow-band) filter (Eq. [5jl and /l e ff,bb(nb) the ef- 
fective wavelength of broad-band (narrow-band) filter (Eq.[4}. 

The effective wavelength of a filter with transmission curve 
T(A) is given by 



teff 



/ AT(A)dA 



j T(A)dA 
and the width of the filter AA by 



(4) 



AA 



I 



T(A)dA/T n 



(5) 



with r max the peak transmission of the filter. For a top-hat filter, 
the effective wavelength is equal to the central wavelength, the 
width equals the FWHM. 

If the central wavelengths of the narrow-band and broad- 
band filters are roughly equal and the Lya line falls in the centre 
of the filters, then eliminating either C or F\jy a by substituting 
Eq. (0 in Eq. Q gives: 



Lya — 



c 



Lya 



A/lbbA/l n b(/^,nb - Abb) 

A/lbb - A/l n b 
A/lbb/.1,bb - A/l n b/,},nb 



A/lbb - A/l n b 

and using Eq. results in an expression for EWohs'. 



obs 



A/lbbA/L n b(/l,nb - Abb) 
A/lbb/i.bb - A/lnb/l.nb 



(6) 
(7) 

(8) 



(e.g. bunker et all 1 19951 iMalhotra & Rhoadsi 12002b . 
Alternatively, Eq. can be used if it is expected that 
the fraction of the continuum flux falling in the filters that is 
absorbed by foreground H i is comparable to the fraction of the 
Lya line that is extinguish ed by intergala ctic absorption (as 
assumed by e.g. lMalhotra & Rhoadsll2002h . 

If the central wavelengths of the narrow-band and broad- 
band filters differ, as is the case with our filters, then the 
slope of the UV continuum is needed to extrapolate the contin- 
uum strength from the central wavelength of the broad-band to 
the central wavelength of the narrow-band. Including an extra 
broad-band contribution redward of the Lya line, the contin- 
uum slope can be calculated as well. Below is described how 
the equivalent width of a z — 3.13 Lya emitter can be computed 
using our available photometry. 

Assume that a Lya emitter has a spectral energy distribu- 
tion that consists of a Lya line with flux f^ya and a UV con- 
tinuum redward of the Lya line with strength C and powerlaw 
slope /3 {f a °c AP). The flux density in the narrow-band (/i, n b), 
V-band (f^y) and /-band (/y) can then be characterized as: 
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Anb - GnbC/i, ff nb + £ n b F]_y a / AA nb 
Av - GvC/lg ff v + eyF^ya/AAy 
Al = C/tf ffI 



(9) 

(10) 
(11) 



with /leff,nb/v/i the effective wavelength corresponding to the 
narrow -band, V and / filter respectively (Eq.0Ji, AA the width of 
the filter (Eq.[5}, e the efficiency of the filter at the wavelength 
of the redshifted Lya line and Q the fraction of the continuum 
flux falling in the filter that is absorbed by the Lya forest (Eq. 
I12l> . It should be mentioned that, in contrast to Eq. ©, no cor- 
rection factor for foreground absorption of the Lya line is ap- 
plied in this calculation. If foreground extinction of the Lya 
line is taken into account, then the equivalent width and Lya 
line flux will be higher by ~ 60% (see Ea. ll3l >. 

For the filters (and instrument) used in this project, the 
input parameters are: /t e ff,nb = 5040.1 A, A^ nb = 61.1 A, 
Aegy = 5561.9 A, AA V = 1145.6 A and A eff j = 7946.5 A. The 
efficiency e of the filters depends on the redshift of the Lya line. 
For all objects a redshift of z — 3.13 is assumed, and the effi- 
ciencies are e nb = 0.76 and ey = 0.74. It should be stressed that 
the computed equivalent width does not depend strongly on the 
assumed redshift in the interval z = 3.13-3.17. Assuming a 
redshift of z = 3.12 will yield equivalent widths that are a fac- 
tor of ~ 2 higher compared to the equivalent widths computed 
wifhz = 3.13. 

The fraction of the continuum flux that is absorbed by fore- 
ground neutral hydrogen averaged over the bandpass is Q: 




Q = 



J e~ TcS T(A)dA 
J T(A)dA 



(12) 



where T e fr is the effective opacity of H i. For observed wave- 
lengths between the redshifted Lya and redshifted Ly/3 line 
(/l Lyj s(l + z) < A obs < A Lya (l + z)), the expression for r eff that q - 
has been taken is: AALQyX 



m v - m nb 

Fig. 3. Color-color diagram for the 3209 objects detected in 
the narrow -band image with a signal-to-noise greater than 5. 
The dashed line shows the color of objects with a rest-frame 
equivalent width of EWo = A. The solid line indicates where 
EWq = 15 A. Objects not detected in the V-band and/or in the 
/-band are plotted with an arrow. 

Subtraction of Eq. (I14> from Eq. dl5> results in an equa- 
tion of the form cP - bP — constant. This equation was solved 
numerically. 

When p was computed, the UV continuum flux density C 
and the Lya line flux F Lya were calculated using Eqs. and 



A^ v - AA' nh f /: 



A,nb 



eff.V 



eff.nb 



r eff = 0.0036 

U216A/ 



3.46 



(13) 



JPress et allll993llMadaulll995l) . Because the Lya line of an 
object at z — 3.13 falls in the blue wing of the V filter, the frac- 
tion of the continuum flux falling in the V filter that is absorbed 
is small and Q is near unity: Qy ~ 0.97. In the narrow-band, 
<2 n b = 0.92 for a source at z = 3.13. 

To calculate the equivalent width of an individual Lya line, 
Eqs. - Q It were solved for /3, C and FL ya . This was done 
in the following way. Equation was multiplied by AA' nb = 
A/inb/^nb and Eq. fllOi by AA'^ = AAy/ey. Substituting C = 
Ai/^ffi (Eq. EJ gave: 



and 



Lya 



Anb/XGnb X 



eff.nb 



) - Av/(GV ^ effV 



l/(AA' nb Q ab A P ^J-l/(AA' v QyA P ee ^ 



(16) 



(17) 



With C and F^ya, the equivalent width (EW) for each object 
was computed: 



C(A hya (l + z)f 



(18) 



A/J| lb Anb - Al AA' nb g„b I eff '" b I + ^Lyu- 



A^v Av = Al AA 'y 2v 



^eff,' 



teff.I 



+ Fi 



Lya- 



(14) 
(15) 



with Ajjy a the wavelength of the Lya line. The rest frame equiv- 
alent width (EWo) is given by: EW = EW ob J(l + z). 

To estimate the uncertainties in the computed parameters, 
the observed flux densities were randomly varied 50 000 times 
over a range having a standard deviation equal to the uncer- 
tainty. The distributions of f3, C, Fi,y a and EWo were used to 
estimate the errors in these quantities. Because the values of 
the equivalent width were not normally-distributed (Gaussian) 
around the central value, two errors were calculated, labelled 
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100 - 




m v - m nb 

Fig. 4. Equivalent width (calculated from Eqs. I9l41 8i versus 
my - 'Wnb color. For clarification, only objects with an EWo > 1 
A are plotted and the error bars are left out. The solid line indi- 
cates a color of my - m n b = 0.72, while the dashed line denotes 
the division between low (EWo < 15 A) and high (EWo > 15 
A) equivalent width objects. 

AEWy and AEW^. These were computed from the values in 
the distribution that were outside the central 99.73% of all val- 
ues. The difference between these values and the central value 
was taken as defining three sigma uncertainties. 

For each object detected in the narrow-band image, the line 
flux, UV continuum and equivalent width and their errors were 
computed. Because no /-band data had been taken yet at the 
time that the candidates for the spectroscopy had to be selected, 
a flat spectrum (J3 = -2) was assumed for all sources. In Fig. 
13 the mi - m-nh color is p lotted against the my - m n \, color. 
Following IVenemans et alJ 12002). objects with EWo > 15 A 
and EWo/AEWy > 3 were selected as candidate Lya emitters. 
Each individual Lya candidate was inspected visually in order 
to remove spurious candidates, like leftover cosmic rays or ob- 
jects in the "spikes" of bright stars. This resulted in a list of 77 
candidate Lya emitters with EWo > 15 A of which 6 had 15 A 
< EWo < 20 A. 

The main difference between the usage of Eq. and Eqs. 
(I9l-jl8t to compute the equivalent width can be seen in Fig. 
@] Using Eq. © (and thereby assuming a fixed slope of /3 = 
-2) the line with EWo = 15 A would lie at a constant my — 
m n b = 0.72 (the solid line in Fig. |4}. As a result, the EWo of 
three very blue objects (with f3 < -2) would be overpredicted, 
falsely selecting objects as Lya emitters (e.g. the three crosses 
in Fig.0}. On the other hand, 1 1 red Lya emitters with /3 > -2 
would not pass the selection criterion my - m n \, > 0.72, while 
their EWo as calculated with Eqs. d9l- (ll8> is greater than 15 A 
(diamonds in Fig.0}. 



spectroscopy session was carried out on 2001, September 22 
with VLT/FORS2 in the multi-object spectroscopy mode with 
19 movable slits of a fixed length of ~ 20 - 22". The night was 
photometric, but because of strong winds (> 12 m s~' ) the see- 
ing fluctuated between 1" and 2". Spectra of 12 candidate Lya 
emitters were obtained for 4 x 2700 s and 1 x 1 800 s through 1 " 
slits with the 1400V grism at a dispersion of 0.5 A pixeL 1 . This 
grism was chosen for a number of reasons. First, it has a high 
peak efficiency of ~85% at wavelengths that corresponds to the 
redshifted Lya line of the radio galaxy. Secondly, because ob- 
servations of high redshift Lya emitting galaxies have shown 
that the width of the Lya line lies predominan t ly in the range 
200-500 kms~' (e.g . IPentericci etafl l2000at bawson et all 
l2002tlHuetalll2004l) . the resolution of the grism (R = 2100, 
corresponding to ~ 150 kms -1 ) ensured that the Lya emis- 
sion line is marginally resolved (see Sect. 15. 1> , maximizing 
the signal-to-noise of the observed line. Also, the resolution is 
large enough to distinguish a high redshift Lya emitting galaxy 
from a low redshift contaminant, the [On] AA3726, 3729 emit- 
ter. With the 1400V grism the [On] doublet is resolved (see 
Fig. [5] for two examples). For the wavelength calibration ex- 
posures of He, HgCd and Ne arc lamps were obtained. The 
spectrophotometric standard st ar LTT 1788 JStone & Baldwinl 
ll983HBaldwin & Stonelll984 was observed with a 5" slit for 
the flux calibration. 

On 2001, October 18, 19 and 20 spectra were obtained with 
FORS2 in the mask multi-object spectroscopy mode. In this 
mode objects are observed through a user defined, laser-cut 
mask with slits which had variable lengths (typically 10 - 12") 
and widths of 1". The nights were photometric with an average 
seeing of 1". The 1400V grism was used to observe 37 candi- 
date Lya emitters in two masks, of which 25 were included in 
both masks. The first mask was observed for 4 x 2700 s and the 
second mask for 10 x 2700 s and 1 x 2100 s. The pixels were 
binned by 2 x 2 to avoid the noise in the spectra being domi- 
nated by read noise. This resulted in a dispersion of 1 ApixeL 1 
and a spatial scale of 074 pixel . Spectra of the standard star 
LTT 1788 were obtained for the flux calibration. 

During the last observing session (2001, November 15 and 
16), the instrument used was FORS1 on Melipal (VLT UT3). 
The main goal of this run was to measure the polarization of 
the radio galaxy (C. De Breuck et al., in preparation). Due to 
constraints on the positioning and orientation of the mask, only 
three candidate Lya emitters could be observed. The width 
of the slits was 1". The total exposure time was 19 800 s. 
The average seeing of these photometric nights was 078. The 
grism used for the observations was the "300V" with a resolu- 
tion of 440, a dispersion of 2.64 ApixeL 1 and a spatial scale 
of 072 pixel -1 . T he spectrophotometric standard stars Feige 
110 a nd LTT 377 l lStone & Baldwinl 1 1983UBaldwin & Stone! 
1984) were observed for the flux calibration. 



4. Spectroscopy 

4.1. Spectroscopic observations 

Spectra of candidate Lya emitters were taken during three sep- 
arate observing sessions (see Table 1 for an overview). The first 



4.2. Data reduction 

The spectra were reduced in the following way. Individual 
frames were flat-fielded using lamp flats, cosmic rays were 
identified and removed and the background was subtracted. 
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4900 4950 5000 5050 5100 5150 

Wavelength (A) 



Fig. 5. Spectra of two [O n] emitters observed in the field (top 
two spectra). For comparison, the spectrum of one of the con- 
firmed Lya emitters is shown at the bottom. The spectra are 
offset from each other by 1.5 x 10~ 18 erg s _1 ctrT 2 A -1 . 



The next step was the extraction of the one-dimensional (ID) 
spectra. Typical aperture sizes were l"-l'/5. If the spectrum 
of the object could be seen in the individual frames, then a 
spectrum was extracted from each frame and these spectra 
were combined. If the object was undetected in the individ- 
ual frames, then the background subtracted two-dimensional 
frames were combined and a ID spectrum was extracted from 
this image. All ID spectra were wavelength calibrated using 
the arc lamp spectra. For spectra taken with the 1400V grism 
the rms of the wavelength calibration was always better than 
0.05 A, which translates to Az = 0.00004 at z ~ 3.13. The 
wavelength calibration with the 300V grism had an rms of 0.8 
A (Az = 0.0007 at z ~ 3.13). A heliocentric correction was ap- 
plied on measured redshifts to correct for the radial velocity of 
the Earth in the direction of the observations. Finally, the spec- 
tra were flux calibrated. The fluxes of the photometric standard 
stars in the individual images were consistent with each other 
to within 5%, so we estimate that the flux calibration of the 
spectra is accurate to ~ 5%. 

4.3. Results 

Spectra were obtained for a total of 40 candidate Lya emit- 
ters, of which 1 1 were observed during two separate observ- 
ing sessions, and the central radio galaxy. Of the 40 candidate 
emitters, only 7 failed to show an emission line. Six of these 
unconfirmed emitters had predicted line fluxes below 10~ 17 
ergs -1 cirT 2 and were probably too faint to be detected. Two 
of the 33 emission line objects showed two lines with almost 
equal strength, separated by ~4 A (Fig.|5}- These objects were 
identified to be [On] /L13726, 3729 emitters at a redshift of 
~ 0.35. One of the [O n] emitters was re-observed in November 
2001 and a nearly flat continuum was revealed with no break 
around the emission line, confirming that the object could not 
be a Lya emitter at z ~ 3.13. None of the other emitters had 
more than one emission line in the spectrum. This excluded 



Table 3. Characteristics of the Voigt absorption profiles derived 
from the spectra. For each absorption profile, its centre relative 
to the peak of the emission line, width (b) and H i column den- 
sity (N) is printed. 



Object Centre b (km s 1 ) log N (cm 2 ) 



344 


-80 = 


b 10 


52 ± 


11 


14.4 ±0.1 


995 


-150 


±20 


74 ± 


59 


16.0 ± 2.5 


1147 


-70 = 


t 20 


101 ± 


24 


14.6 ±0.1 


1203 


-90 = 


t 10 


104 ± 


11 


14.9 ±0.1 


1518 


80 ± 


50 


38 ± 


16 


13.8 ±0.2 




-210 


+ 60 


155 ± 


52 


15.0 ±0.2 


1612 


-150 


± 10 


79 ± 


16 


14.4 ±0.1 


1710 


-130 


± 20 


80 ± 


30 


14.1 ± 0.2 


1867 


-60 = 


t 20 


30 ± 


16 


13.1 ±0.3 


2487 


250 ± 


170 


108 =t 


14 


14.6 ±0.1 




-1150 


±200 


629 ± 


93 


16.0 ± 0.2 


3101 


-40 ± 


110 


193 ± 


131 


14.4 ± 0.4 




-240 


±20 


79 ± 


38 


14.3 ± 0.3 


3388 


-130 


± 10 


59 ± 


13 


14.3 ±0.1 


HzRG 


200=1 


= 10 


151 : 


t9 


14.9 ±0.1 




-270 


± 10 


245 ± 


25 


14.9 ±0.1 




-660 


± 10 


144 ± 


20 


14.8 ±0.1 




-970 


±20 


131 ± 


27 


14.2 ±0.1 



identification of the emission line with [Om] /15007, because 
then the confirming [Oiii] /14959 would have been visible. 
Furthermore, a number of emitters showed an asymmetric line 
profile (Figs. ll5H20t . a feature often seen in spectra of high red - 
shift Lya emitters (e.g. lAiiki et"aHl2002tlDawson et alll2002l) . 
Therefore, the 3 1 remaining emitters were identified with be- 
ing Lya emitters. The fraction of contaminants in our sample is 
2/33 = 6.1%, similar to the fraction of low redshift interlopers 
of 6.5% in the study of Lya emitters at z ~ 3.09 of St eidel et all 
J200(j) . 

5. Properties of the Lya emitting galaxies 

The one-dimensional Lya emission lines were fitted by a 
Gaussian function and - if absorption was clearly present - 
in combination with Voigt absorption profiles. The best fit 
Gaussian was used to calculate the redshift, line flux and 
FWHM of each emitter. In Table 2 the properties of the con- 
firmed Lya emitters are summarized. The IDs correspond to 
the object's number in the SExtractor catalog. The rest-frame 
equivalent width EWo was taken from the imaging. The star 
formation rate (SFR) was calculated using the Lya line flux de- 
rived from the images, and assuming Case B recombination and 
using the Ha luminosity to SFR conversion from M adau et alJ 
1 19981 see Sect.ES. 

5.1. Line profiles 

As mentioned in the previous paragraph, emitters which clearly 
showed an emission line with an absorbed blue wing (see Figs. 
I15H20> were fitted by a Gaussian emission line with one or 
more Voigt absorption profiles. The characteristics of the ab- 
sorption profiles are listed in Table 3. The other emitters were 
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Table 2. Position and properties of the Lya emission line of the 3 1 confirmed Lya emitters and the radio galaxy. 
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2, 

J 


1 A07 




0007 
U.UUU / 


1.7 ± 


0.4 


TT + 12 


600 + 


130 


n q + 


1 


1968 


03 


18 05.14 


-25 34 51.2 


3 


1564 


+ 


0.0004 


0.4 + 


0.1 


< 66 


300 + 


80 


1.7 ± 


0.3 


2413 


03 


18 20.40 


-25 33 00.3 


3 


1339 


+ 


0.0002 


0.7 + 


0.2 


62+?? 

-11 


200 + 


50 


2.7 ± 


0.3 


2487 


03 


17 59.63 


-25 34 01.3 


3 


1175 


± 


0.0024 


19.3 d 


:2.6 


202+' 


2510 + 


190 






2637 


03 


18 22.52 


-25 34 10.1 


3 


1224 


± 


0.0002 


1.3 ± 


0.3 


7Q+43 
/y -l6 


200 + 


50 


3.2 ± 


0.5 


2871 


03 


18 25.71 


-25 33 40.5 


3 


1334 


± 


0.0003 


0.4 + 


0.1 


>75 


250 + 


70 


0.7 ± 


0.1 


3101 


03 


18 24.09 


-25 32 12.2 


3 


1313 


± 


0.0003 


4.3 + 


1.2 


114+47 


800 + 


100 


5.2 ± 


0.3 


3388 


03 


18 26.09 


-25 32 53.1 


3 


1379 


± 


0.0002 


2.1 ± 


0.1 


112+ 26 


370 + 


20 


8.6 ± 


0.3 


HzRG 


03 18 12.01 


-25 35 10.8 


3.1307 


± 


0.0001 


155.1 


±4.2 


257- 


1320=! 


: 10 


76.1 ± 


0.4" 



" Object contains an AGN, therefore the estimated SFR is unreliable. 

fitted by a single Gaussian. The observed width of the line 
was deconvolved with the instrumental width, which was 150 
kms -1 (see Sect. 14. 1>. 

The radio galaxy has an emission line that can be fitted by 
a Gaussian with a FWHM of ~ 1300 kms' 1 . The line w idth 
is very similar to that of other HzRGs (e.g. iDe Breuck et all 
l200ll IWillott et all 120021) . Only one of the confirmed emit- 
ters has a broad Lya line. Emitter #2487 has a line FWHM of 
~2500 kms -1 , and is therefore likely to also harbour an AGN. 
The FWHM of the Lya emission line of the rest of the emitters 
ranges from 120 kms -1 to 800 kms"' (Fig.|6jl, with a median 
of 260 km s~' and a mean of 340 km s _1 . 

The inferred column densities of the absorbers are in the 
range of 10 13 — 10 16 cm" 2 (see Sect.|6}. Using the spatial extend 
in the 2D spectra as an estimate of the size of the H 1 absorber, 
the amount of projected neutral HI near the emitters is in the 
range of > 2 x 10 2 - 5 x 10 4 Mq (see Sect.|6]for the details of 
the individual emitters). For the fainter Lya- emitters, it cannot 
be excluded that the troughs are due to substructure in the Lya- 
emitting regions, rather than H 1 absorption. 
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Fig. 6. Histogram of the line widths of the emitters. 
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Fig. 7. / magnitude versus UV continuum slope [3. The dashed 
line corresponds to/? = -2.1, the color of an unobscured con- 
tinuously star forming galaxy with an age of ~ 10 s — 10 9 yr. 
The hashed area indicates the color of a young (~ 10 6 yr), star 
forming galaxy. 



5.2. Continuum colors 

In Table 4, the magnitudes and UV continuum slopes of the 
confirmed emitters are listed (for a description on how the slope 
was calculated, see Sect. 13. 31 . For objects detected in the ACS 
image which were not resolved into several components, the 
Ig 14 magnitude and the continuum slope /?acs calculated using 
this magnitude is given. In Fig. Q the / magnitude is plotted 
against the continuum slope. Excluding the radio galaxy and 
emitter #2487, which contains an AGN (Sect.l5~Tt. the UV con- 
tinuum slope of the confirmed emitters ranges from f3 = 0.62 
to/3 = -4.88 with a median of /3 = -1.76. 

The blue median color of the Lye emitters may be due to a 
selection effect. The candidate emitters for spectroscopy were 
selected when only one broad-band flux was available and a 
slope of /3 = -2 was assumed for all objects to compute the 
equivalent width (see Sect. l3.3l and Fig. [4}- Because the narrow- 
band is on the blue side of the broad-band filter that was used, 
the equivalent width and line flux of bluer objects with f3 < 
-2 tend to have been overestimated, while "red" objects (J3 > 
-2) have an equivalent width and Lya flux that are likely to be 
underestimated (see Sect. I3.3I >. This effect could have biased 
the spectroscopic sample towards blue objects. For example, 
the emitter with the bluest color, #1446, has an equivalent width 
of EW{) — 12 A, which falls below the selection criteria (EWq > 
15 A), but the object was selected for spectroscopy because 
it had ;«y - m n \, > 0.72 (see discussion at the end of Sect. 
13.31 . To determine the effect of this bias, the color of a flux 
limited sample was determined. There are 3 1 (candidate) Lya 
emitting objects with a Lyor flux > 1.5 x 10" 17 ergs -1 cm"" 2 
in the field, of which 25 are confirmed. Again excluding the 
radio galaxy and the AGN, the median color of the remaining 
29 emitters is /3 — - 1 .70. This is bluer than the average color 
of Lya- emitting Lyman Break Galaxies ( LBGs) , which have a 
slope of/3 = -1.09 + 0.05 (Ehapley et all l2003l) . 



Models of galaxies with active star formation predict UV 
continuum slopes in the range f3 = -2.6 to f3 — -2.1 for an un- 
obscured, continuously star formi ng galaxy with ages between 
a few Myr and more than a Gyr llLeitherer etal l ll999l) . 18 out 
of the 27 (67%) confirmed Lya- emitters for which f3 could be 
measured, have colors within 1 cr consistent with being an un- 
obscured starburst galaxy. Of those, 15 (56% of the sample) 
have such blue colors with 1 cr that they could be star forming 
galaxies with ages of order 10 6 yr, which have f3 2.5. 



5.3. Morphologies 

Of the 32 confirmed Lya- emitting sources, 19 (including the 
radio galaxy) were located in the area that was imaged by the 
ACS (Fig. I14> . Two of these emitters remained undetected to 
a depth of Igu > 27.1 mag arcsec~ 2 (3 <x). On the position of 
the radio galaxy, the ACS image shows several objects within 
3" (~25 kpc), surrounded by low surface brightness emission 
(> 24.8 mag arcsec~ 2 , see Fig. [8). Such a clumpy structure 
is oft en seen at the position of HzRGs (e.g. IPentericci et all 
Il999l) . Interestingly, there are three other Lyor emitters with 
morphologies that resemble the radio galaxy (Fig. [HJ. Each of 
these three objects consists of at least three clumps of emission, 
which are less than one kpc separated from each other. The re- 
mainder of the confirmed emitters can be identified with single 
objects in the ACS image. 

To quantify the size of these objects, the half light ra- 
dius (ri,) of each emitter was measured using the program 
SExtractor. The half light radius is defined as the radius of 
a circular aperture in which the flux is 50% of the total flux. 
However, as already discussed in Sect. 13.21 and shown in Fig. 
12 the fraction of the total flux of an object that is missed by 
SExtractor increases when the object is fainter. As a conse- 
quence, the half light radius that is measured by SExtractor 
would underestimate the size of the object. To determine how 
much the half light radius was underestimated, galaxies with 
a range of sizes were varied in brightness and added to the 
ACS image, and the half light radii of those objects was mea- 
sured by SExtractor. It was found, as mentioned above, that 
the fainter the object, the smaller its measured r;,, an effect that 
was stronger for larger galaxies, see Fig. ^2 Using the results 
of these simulations, an attempt could be made to correct the 
measured sizes of the confirmed Lya emitters. Unfortunately, 
this correction could overestimate the true size of compact ob- 
jects (i.e. objects with a half light radius similar to that of stars). 
However, this only strengthens our conclusions (see below). In 
Table 5 the sizes of the emitters in the /814-band are printed. 
The half light radii of the emitters range from 0'.'06 to 0V18. 
The error in the half light radius is defined as the half light ra- 
dius divided by the signal-to-noise of the object. Translating 
the sizes directly to physical sizes, the measured half light radii 
correspond to 0.5-1.5 kpc. The median size is ~1 kpc. 

The mean half light radius of isolated, unsaturated stars in 
the field was found to be ~0706. Four of the emitters in the 
ACS field have a r/, that is within 1 cr equal to the half light 
radii of the stars in the field. These four emitters are classified 
as unresolved (Fig.|5}. 
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Table 4. Narrow-band, V, I and ACS 7gi4 magnitudes and UV continuum slopes ft of the confirmed emitters. Colors were 
measured in a circular aperture, while elliptical apertures were used to determine the total magnitudes (see Sect. 13. 2\ . For objects 
with a signal-to-noise less than two a 2 cr upper limit is given. 



Object NB magnitude V magnitude / magnitude /gi4 magnitude /3 /3acs 



344 


24.63 ± 


0.06 


26.62 ± 


0.22 


26.78 ± 0.28 


27.76 ± 


0.23 


-1.43 ± 1.21 


-3.91 ± 


1.06 


695 


24.89 ± 


0.09 


26.42 ± 


0.20 


25.96 ± 0.19 


26.53 ± 


0.10 


0.62 ± 1.05 


-0.92 ± 


1.09 


748 


25.45 ± 


0.14 


27.03 ± 


0.37 


27.13 ± 0.42 


— 




-1.75 ± 1.85 


— 




995 


24.97 ± 


0.08 


26.14 ± 


0.16 


26.49 ± 0.23 


26.59 ± 


0.11 


-2.40 ± 0.72 


-2.65 ± 


0.63 


1029 


25.00 ± 


0.07 


26.19 ± 


0.12 


26.76 ± 0.23 


— 




-3.18 ± 0.74 






1099 


25.23 ± 


0.13 


26.39 ± 


0.28 


26.49 ± 0.28 


— 




-1.98 ± 0.96 






1147 


25.74 ± 


0.12 


27.83 ± 


0.44 


> 28.38 






< -1.40 






1203 


24.79 ± 


0.06 


26.00 ± 


0.11 




27.30 ± 


0.15 


— 


-4.77 ± 


0.55 


1361 


24.76 ± 


0.10 


26.47 ± 


0.24 


26.63 ± 0.27 


— 




-1.65 ± 1.14 






1395 


26.40 ± 


0.21 


27.50 ± 


0.47 


> 27.69 






< -2.30 






1446 


25.53 ± 


0.10 


26.32 ± 


0.12 


27.46 ± 0.29 


26.55 ± 


0.12 


-4.88 ± 0.96 


-2.46 ± 


0.46 


1498 


26.14 ± 


0.14 


27.55 ± 


0.29 


27.79 ± 0.50 


27.12 ± 


0.18 


-2.20 ± 2.09 


-0.49 ± 


1.06 


1518 


23.65 ± 


0.02 


24.47 ± 


0.03 


24.54 ± 0.04 


24.57 ± 


0.03 


-2.20 ±0.11 


-2.12 ± 


0.12 


1612 


25.79 ± 


0.16 


27.65 ± 


0.54 


> 27.38 


27.84 ± 


0.28 


— 


> -0.65 


1710 


25.51 ± 


0.08 


27.09 ± 


0.19 


27.37 ± 0.38 


27.34 ± 


0.19 


-2.26 ± 1.40 


-2.18 ± 


0.78 


1724 


25.74 ± 


0.15 


> 27.91 


> 27.44 


~~ ^ 




— 






1753 


24.72 ± 


0.09 


25.67 ± 


0.18 


25.71 ± 0.20 


— C 




-1.90 ± 0.60 






1759 


24.77 ± 


0.09 


25.87 ± 


0.19 


26.02 ± 0.27 


— 




-2.03 ± 0.83 






1829 


25.11 ± 


0.13 


25.73 ± 


0.26 


25.02 ± 0.25 






-0.07 ± 0.60 






1867 


22.86 ± 


0.01 


24.25 ± 


0.03 


24.17 ± 0.04 


— ° 




-1.42 ± 0.15 






1891 


25.56 ± 


0.13 


26.76 ± 


0.29 


> 27.24 


27.45 ± 


0.20 


< -3.02 


-3.50 ± 


1.01 


1946 


25.58 ± 


0.18 


26.40 ± 


0.36 


> 26.76 






< -2.93 








Ij.jy ± 


V. 1 J 


27.56 ± 


0.47 


Z/.Zj ± U.j4 


Z/.jy) ± 


r\ 1 o 
U. to 


> — U.jz 


> -1.03 


1962 


25.44 ± 


0.08 


26.44 ± 


0.11 


26.07 ±0.13 


26.10 + 


0.08 


-0.85 ± 0.48 


-0.99 ± 


0.38 


1968 


25.00 ± 


0.14 


26.48 ± 


0.45 


> 26.37 


26.61 ± 


0.12 


< -1.23 


-1.89 ± 


1.63 


2413 


24.51 ± 


0.07 


25.97 ± 


0.16 


26.04 ± 0.23 






-1.76 ± 0.73 






2487 


22.25 ± 


0.01 


23.40 ± 


0.02 


23.30 ± 0.02 






-0.57 ± 0.09 






2637 


24.88 ± 


0.09 


25.97 ± 


0.17 


25.88 ± 0.20 






-1.23 ± 0.64 






2871 


26.04 ± 


0.18 


> 27.64 


> 27.22 












3101 


23.99 ± 


0.04 


25.70 ± 


0.12 


25.49 ±0.15 






-0.71 ± 0.60 






3388 


23.31 ± 


0.03 


25.25 ± 


0.09 


25.54 ± 0.16 






-1.92 ± 0.52 






HzRG 


21.19'' ± 


: 0.01 


23.31'' ± 


:0.02 


22.97 rf ± 0.02 


c 




0.27 ± 0.09 







" Photometry unreliable due to nearby bright star. 
* Undetected in the ACS image, / 814 > 27.1 mag arcsec 
c Resolved by HST into several components. 
d Photometry influenced by nearby objects. 



The sizes of the confirmed Lya emitters can be compared 
to other high redshift galaxies, e.g. LBGs. Recently, sizes 
were measured of galaxies at various redshift s in the Great 
Observatories Origins Deep Survey (GOODS, Ferguson et al., 
12004 . For their analysis, they used SExtractor with circular 
apertures having a radius that is 10 times larger than the first 
radial moment of the light distribution to ensure that all the flux 
was inside the aperture. The survey was restricted to rest-frame 
luminosities between 0.7 L* and 5 L*. Us i ng the luminosity 
function of LBGs derived by ISteidel et alJ dl999h . this corre- 
sponds to a magnitude range of 22.78 < m.R < 24.92. Only two 
confirmed emitters located in the ACS field (#1518 and #1867) 
satisfy the luminosity criterion used in the GOODS analysis. 
Emitter #1867 is resolved into several clumps of emission. The 
size of emitt er #1518 meas ured with the same input parame- 
ter aslFer guson et alJll2004 . is 07106 ± 07006, consistent with 



the 07102 + 07003 derived using our own input parameters. 
The half light radius of emitter #1518 is among the smallest 
Ferguson et al. are finding. The average size of LBGs at z ~ 3 
is 0728 (~ 2.3 kpc). Thus, the Lya emitters are small compared 
to LBGs at the same redshift, provided that the method we used 
to measure the sizes of the Lya emitting galaxies gives compa- 
rable half light radii as the approach of Ferguson et al. (as was 
the case for emitter # 1518). 

5.4. Star formation rate 

The average star formation rate (SFR) of the confirmed emit- 
ters, as derived from the Lya flux (see Table 2), is 2.5 Mq yr~' 
(excluding the radio galaxy and the QSO, emitter #2487). This 
calculati on assumed a Lyq/H a ratio of 8.7 (Case B recom- 
bination, Brocklehurst, Il971l) and a Ha luminosity to SFR 
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Table 5. Half light radii of the confirmed emitters located 
within the field of the ACS. 



Object 


n (") 


r h (kpc) 


s/n 


344 


0707 ± 0703 


<0.8 


4.8 


695 


0710 ±0702 


0.6 ± 0.2 


10.5 


995 


a.' 09 ± 0701 


0.5 ± 0.2 


10.2 


1203 


a.' 06 ± 0701 


<0.6 


7.1 


1446 


0714 ± a' 02 


1.0 + 0.2 


9.4 


1498 


0718 ±0706 


1.4 ±0.5 


6.0 


1518 


a' io ± a'oi 


0.7 ±0.1 


43.1 


1612 


0714 ±0705 


1.0 ±0.5 


3.8 


1710 


o.'io±a'03 


0.6 ± 0.4 


5.8 


1753 


Clumpy 


— 


~ 12 


1829 


Clumpy 




~ 7 


1867 


Clumpy 




~ 20 


1891 


0710 ±0703 


< 1.3 


5.5 


1955 


0708 ± ff.'03 


< 1.0 


6.1 


1962 


0.' 13 ± 0701 


1.0 ±0.1 


12.8 


1968 


0714 ± a' 02 


1.1 ± 0.2 


9.0 


HzRG 


Clumpy 




~ 35 



conversion for a S alpeter initial mass function (IMF) from 
lMadauetalJ (fr998): 



SFR 



-'Ha 



Ha 



1.6 x 10 41 ergs- 1 



(19) 



Because of Lya absorption (see e.g. Fig. ll5H20l . this SFR cal- 
culation gives a lower limit. 

An alternative way to estimate the SFR is to use the level of 
the UV continuum. The flux density at a wavelength of A KSt = 
1500 A can be converted to a SFR following the relation 



S FR\jy - 



AjvOW = 1500 A) 
8.0 xlO 27 ergs- 1 Hz" 1 



(20) 



for a Salpeter IMF jMadau et all Il998l) . In Fig. El the SFR 
as calculated from the Lyor emission is plotted against the UV 
continuum SFR. On average, the two methods to calculate the 
SFR give the same result, with the SFR measured from the UV 
continuum a factor 1-1.5 higher than the Lya inferred SFR. 
The average SFR of the emitters as measured by the UV con- 
tinuum is < 3.8 M©yr -1 . This is much lower than the aver- 
age SFR of LBGs, whic h is somewhere between 10 and 100 
M yr-' (e.g. lGiavaliscoll2002l) . 

Recent measurements of the polarization of the UV con- 
tinuum of the radio galaxy indicate that the UV continuum 
is dominated by emission from stars. The contribution from a 
scattered AGN is small, which is implied by the upper limit on 
the polarization of the continuum of P < 4% (C. De Breuck et 
al., in preparation). If all the light at a rest-frame wavelength of 
1500 A is due to young stars, then the SFR of the radio galaxy 
is 40.5 + 0.8 Mq yr . No correction is made for dust absorp- 
tion. This is similar to the uncorrected SFR (as calculated from 
the rest-frame UV continuum) in radio galaxies at z ~ 2.5 (e.g. 
IVernet et all 1200 lh and a factor of ~ 5 lower than the SFR of 
the radio galaxy 4C41.17 at z = 3.8 JChambers etail Il990t 
iDev et allll997h . 



6. Notes on individual objects 



- #344: The spectrum of this emitter is shown in Fig.^] The 
Lya line can be fitted by a combination of a Gaussian and 
a Voigt absorption profile which is located 80 ± 10 kms -1 
blueward of the peak of the Gaussian. The fit is shown as 
the solid line in Fig.^] From the 2D spectrum a lower limit 
of 2" (~ 16 kpc) on the linear size of the absorber could be 
derived, giving a lower limit on the H i mass of M(H i) > 
550 M Q . 

- #995: The spectrum of this object can be fitted by a 
Gaussian with a narrow Voigt profile centred 150 + 20 
kms 1 blueward of the emission peak (solid line in Fig. 
I15> . The column density of the absorbing neutral hydrogen 
is nearly unconstraint by our spectrum and lies in the range 
JV(Hi) ~ 10° 518 5 . The absorber has a size of > 2" in the 
2D spectrum and taking a column density of 10 16 cm 2 , the 
inferred mass of neutral hydrogen is at least 2 x 10 4 Mq. 

- #1147: This object is one of the two emitters that were 
not detected in the ACS image down to a 3 sigma mag- 
nitude limit of 27.1 magarcsec -2 . Its redshift is ~ 2400 
kms 1 redder than the median redshift of the emitters, and 
is therefore unlikely to be associated with the protocluster 
(see Sect. l7.2> . The spectrum shows absorption which is lo- 
cated 70±20 km s _1 blueward of the peak of the unabsorbed 
emission (Fig. [16}. A lower limit of 550 Mq can be given 
for the H i mass. 

- #1203: This emitter is unresolved in both the VLT and ACS 
images (Fig.[9j. The spectrum is shown in Fig. [HO The ab- 
sorption is located 90 ± 10 kms 1 to the blue of the peak 
of the Gaussian. We obtain a lower limit for the mass of 
neutral hydrogen responsible for the absorption of > 1700 
Mq. 

- #1446: The emitter has a very blue continuum slope (J3 = 
-4.88+0.96). The computed equivalent width of EW = 12 
A is below the selection criterion (EWo > 15 A). The VLT 
/-band magnitude and the magnitude derived from the 
ACS image are different at the 2.8 <x level, taking the dif- 
ferences in effective wavelength of the filters into account. 
Using the ACS magnitude the continuum slope becomes 
-2.61 and the EWo ~ 20 A. The object is resolved in the 
ACS image and has a half light radius of 1 .0 + 0.2 kpc (Fig. 
EH). 

- #1498: In the ACS image, an object with r h = 1 .4 + 0.5 kpc 
is located ~075 from the position of the weak emitter in the 
VLT image. 

- #1518: This is the fourth brightest Lya- emitter in the field. 
The object is extended in both the Lya image and the ACS 
image (see Sect. l5.3l and Fig. I10i. The emission line can be 
fitted by a Gaussian with two Voigt profiles superimposed, 
one 80 ± 50 km s -1 to the red and the other 210 + 60 km s 1 
to the blue of the Gaussian (Fig.ll7>. This results in a lower 
limit of the mass of H i of 2300 Mq. 

- #1612: This emitter has a faint continuum (7si4 = 27.84 ± 
0.28) and is barely detected in ACS image (signal-to-noise 
of ~4), and is marginally resolved (Table 5, Fig. II 0i . The 
spectrum shows absorption of io I4 - 4±ai cirr 2 Hi, located 
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150 ± 10 kms 1 to the blue of the redshift of this galaxy 
(FigJQ. 

- #1710: This is a blue emitter (J3 = -2.26 ± 1.40) with an 
absorption trough on the blue wing of the emission line (see 
Fig. ITS), the result of at least 200 Mq of H i. 

- #1753: At the position of this Lya emitter, the ACS image 
shows three separate objects located within ~ 8 kpc (Fig. 

E 

- #1829: This object is resolved by the ACS into an elongated 
structure consisting of several objects (Fig. ID' 

- #1867: Denoted as galaxy "A" by LF96, a spectrum of this 
object taken under bad seeing conditions confirms the red- 
shift measured by LF96 (Table 2). The Lya line is asym- 
metric and can be fitted by a Gaussian and one Voigt ab- 
sorber, which is 60 + 20 kms -1 blueward of the Lya peak 
(Fig.ll8>. The VLT narrow-band image shows an extended 
Lya halo of ~25 kpc (Fig.|D, while in ACS image the ob- 
ject is very clumpy. 

- #1968: This emitter was undetected in the VLT /-band, but 
in the ACS image an object with a half light radius of 1 . 1 ± 
0.2 kpc is visible. 

- #2487: This emitter has the brightest Lya line in the field 
after the radio galaxy, and is called galaxy "B" in LF96. 
As mentioned in Sect. 15. II the Lya line is broad (FWHM ~ 
2500 km s -1 ) which is most likely caused by an AGN. The 
spectrum is characterized by a large absorption trough with 
a column density of 7V(Hi)~ 10 14 6 cm" 2 . Furthermore, the 
red wing of the Lya line is much broader than the blue 
wing. The spectrum can be fitted with two absorbers, lo- 
cated 250 ± 170 kms" 1 to the red and 1150 ± 200 kms" 1 
to the blue of the centre of the emission line. The inferred 
mass of Hi is > 5 X 10 4 Mq. 

- #2719: This is galaxy "C" from LF96. Galaxy "C" was not 
selected as a candidate Lya emitter in our images. It has 
colors comparable to those quoted in LF96, but an EWq of 
l.O^J'j. LF96 measured an EWo > 12 A and a line flux of 
~ 5 x 10 -17 ergs -1 cm -2 , although no spectrum was taken 
of this object to confirm the existence of the line. An ex- 
planation for the fact that this galaxy is not selected by us 
as an emission line candidate could be that the large width 
of the narrow-band filter used by LF96, making it sensi- 
tive to a wider redshift range than our filter. Their filter 
was sensitive to Lya emitters having redshifts in the range 
Z = 3.08 - 3.16, while our filter is sensitive to the redshift 
range z = 3.12 - 3.17. Galaxy "C" could be a Lya emitter 
with a redshift between z — 3 .08 and z = 3 . 12, and it would 
therefore be part of the protocluster, but not be included as 
one of our candidates (see Sect. 17.21 . 

- #3101: The Lya line of this emitter is broad (800 ± 100 
km s -1 FWHM, see Fig.ll9>. as compared to the median line 
width of the emitters (260 kms -1 ). The spectrum can be 
fitted by a Gaussian, superimposed by two Voigt absorbers 
located 40 ± 110 and 240 ± 20 kms -1 to the blue of the 
emission, the result of at least ~ 1000 Mq of H i. 

- #3388: This blue Lya emitter (/? = -1.92 ± 0.52) shows an 
absorption trough 130+ 10 km s -1 blueward of the emission 
redshift (Fig.i20L implying a neutral Hydrogen mass of > 
625 M Q . 



- Radio galaxy MRC 0316-257: The absorption structure of 
the radio galaxy is complicated. Only a Gaussian emission 
line profile with 4 separate absorbers gives a reasonable fit 
(solid line in Fig.[2U|and Table 3). The absorbers are 200 + 
10 km s -1 to the red of the peak of the Gaussian and 270± 10 
kms -1 , 660 ± 10 kms -1 and 970 + 20 kms -1 to the blue. 
Approximately 2" to the north-east of the radio galaxy is a 
foreground galaxy, with a clear spiral structure in the ACS 
image. An emission line was detected in a spectrum of this 
object, with a wavelength around 6965 A, most likely [O n] 
at z ~ 0.87 (C. De Breuck, private communications). 

7. A protocluster at z = 3.13? 
7.1. Volume density 

Several surveys have been carried out to estim ate the (field) 
volume density of Lya emitters at z ~ 3 (e .g. |C^wje_&Hul 
1998 1 iHu et all Il998t iKudritzki et all l200Cj ICiardullo et all 
2002: lFvnbo et alll2003l) . To estimate the (over)density of Lya 
emitters in our field, we will compare our n umbers with those 
found in the survey of ICiardullo et all (I2002[). since it cover s 
the largest area of all the surveys, and of iFvnbo et all (J2003), 
because they present the deepest and spectroscopically most 
com plete comparison sa mple of blank field z ~ 3 Lya emitters. 

ICiardullo et aTl £2002) made a blank field survey to esti- 
mate the density of emission line sources and to calculate the 
contamination of intra-cluster planetary nebulae searches. They 
searched for faint emission line sources in a 0.13 deg 2 field at 
a wavelength of 5019 A. They found 21 objects with an ob- 
served equivalent width greater than 82 A and a m„b < 24.3. 
Assuming all their sources are Lya emitters at z ~ 3.13, the vol- 
ume density of field emitters is = 2.4+9j x 10~ 4 Mpc~ 3 . 
Applying the same selection criteria to our data, 5 emission 
line objects (excluding the radio galaxy) are found in the field 
around 0316-257, all confirmed to be Lya emitters at z ~ 3.13. 
This gives a density of «03i6 = 5.4*H x 10~ 4 Mpc~ 3 . The den- 
sity in the 03 16 field is therefore a factor 2. 2*Jq higher than the 
field density. The large errors on this number are due to small 
number statistics. 

More recently, IFvnbo et all d2003h presented the first re- 
sults of a program to detect faint Lya emitters at z ~ 3. They 
used the same VLT narrow-band filter as described in Sect. 12. II 
to image a field that contained a damped Lya absorber. The 5cr 
detection limit for point sources in their narrow-band image is 
m„b < 26.5 as measured in a circular aperture with a size twice 
the seeing FWHM, is very similar to our 5<x detection limit 
(m„b < 26.4). They found 27 candidate Lya emitters with an 
equivalent width greater than 12.5-25 A, the limit depending 
on the predicted line flux. Subsequent spectroscopy confirmed 
that 18 of the 22 candidate emitters observed are Lya emit- 
ters and two were foreground [On] emitters. Assuming that 
the seven unconfirmed candidate emitters are all Lya emitters 
at z ~ 3, the number of Lya emitters down to a flux limit of 
6 X 10~ 18 erg s -1 cm -2 in their field is 25. Our number of emit- 
ters selected with the same equivalent width limits is ~75 after 
correction for foreground contaminants. This implies a density 
of 3.0+y 7 times the field density. Roughly, we find three times 
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the number of Lya emitters as might be expected from field 
sur veys. 

iMaier etail d2003l) gathered measured abundances of Lya 
emitters from the literature, shifted them to z = 3.5 and fitted a 
model function through the points (a description o f the model 
can be found in lThommes & Meisenheimenl2004|) . They pre- 
dict approximately 2325 Lya emitters per deg 2 in a volume 
with Az = 0.1 with line fluxes exceeding 5 xlO -1 8 ergs -1 cm -2 . 
If their model is correct, then ~ 15 Lya emitting galaxies at 
Z = 3 . 13 are expected within our volume brighter than 7 x 1 0~ 1 8 
ergs -1 cm -2 . Applying this limit, we find 63 galaxies or 59 
galaxies if we correct for possible foreground contaminants. 
Of these, 29 are spectroscopically confirmed. The density is a 
factor 4.0^J 5 higher than the model prediction, in agreement 
with the above estimates. 

To summarize, the density of Lya emitters near the radio 
galaxy is a factor 2-4 higher than the field density, indicating 
that the radio galaxy might reside in an overdense region. 

7.2. Velocity distribution 

The velocity distribution of the emitters is plotted in Fig. [O] 
The response of the narrow-band filter used to select the candi- 
date emitters for spectroscopy is also shown. Interestingly, the 
emitters are not distributed homogeneously over the filter, but 
most of them appear to cluster on the blue side of the filter. The 
emitters which show absorption in their emission line profiles 
seem to be distributed more homogeneously, but this could be 
due to the small number of objects. 

To test whether the clustering of emitters in redshift space 
is significant, Monte Carlo simulations of the redshift distribu- 
tion were performed. 10000 realizations with 31 emitters were 
reproduced, with the narrow-band filter curve as the redshift 
probability function for each emitter. The mean of the observed 
redshift distribution (z = 3.136) differs 2.6 cr compared to the 
simulated distribution (z = 3.146 + 0.004) and the width of the 
observed redshift distribution differs by 1 .7 cr (0.012 compared 
to 0.022 + 0.006). In total, the measured redshift distribution 
deviates from the simulated one by 3.07 cr. This means that a 
redshift distribution as observed was reproduced in only 0.2% 
of the cases. The peak of velocity distribution of the Lya emit- 
ters lies within 200 kms -1 of the redshift of the radio galaxy 
(Fig. ll 3I >. providing evidence that most of the Lya emitters are 
physically associated with the radio galaxy. Taking together, 
the observed overdensity of Lya emitters in our field (Sect. 
17. 1> . combined with the peak in the redshift distribution pro- 
vide compelling evidence that the Lya emitters reside in a pro- 
tocluster atz ~ 3.13. 

8. Properties of the protocluster 

In this Section, we discuss the macro properties of the Lya 
emitters in the protocluster. 



8.1. Velocity dispersion 

To determine the velocity dispers ion of the e mitters, the bi- 
weight scale estimator was used jBeers et all Il990l) . This is 



the most appropriate scale estimator for samples of 20-50 ob- 
jects (Be ers et all ll990). The velocity dispersion is 640 ± 195 
kms -1 , corresponding to a FWHM of 1510 ± 460 kms -1 . This 
is significantly smaller than the width of the narrow-band fil- 
ter, which has a FWHM of ~ 3500 km s -1 . Although most Lya 
emitters are likely to be members of the protocluster, the three 
emitters with velocities > 1500 km s -1 from the peak of the dis- 
tribution (Fig. 1131 are probably field galaxies. Ignoring these 
three field galaxies, the velocity dispersion drops to 535 ± 100 
km s -1 . On the lower redshift side (negative velocities), no clear 
edge is visible in the distribution. This could be due to the de- 
crease in sensitivity of the narrow-band filter on this side of the 
redshift distribution. If the protocluster extends to much lower 
redshifts, our estimate of the velocity dispersion is a lower 
limit. 

8.2. Spatial distribution 

The spatial distribution of the emitters is shown in Fig. El 
where the circles, diamonds and triangles represent the spec- 
troscopically confirmed emitters, and with the sizes of the sym- 
bols scaled according to the velocity offset from the median of 
the emitters. The pluses are unconfirmed candidate Lya emit- 
ters satisfying our selection criteria (see Sect. 13. 3> . The major- 
ity of these candidates (96%) has not yet been observed spec- 
troscopically, while the remaining 4% were too faint to be con- 
firmed. The imaging field of view (3.3 x 3.3 Mpc 2 at z — 3.13) 
is not large enough to show clear boundaries of the structure. 

8.3. Mass 

At a redshift of z = 3.13, the age of the Universe is only 2.2 
Gyr. Taking the velocity dispersion as a typical velocity for a 
galaxy in the protocluster, it would take at least 5 Gyr to cross 
the structure, making it highly unlikely that the protocluster is 
near virialization. Therefore, the virial theorem cannot be used 
to calculate the mass of the protocluster. 

Another way to compute the mass is to use the (comov- 
ing) volume V occupied by the overdensity, the (current) mean 
density of the Universe p and the mass overdensity of the pro- 
tocluster <S„, : 



M = pV(l+S m )=pV{l + S gai /b). 



(21) 



where b is the bias parameter (b = 6 ga i/5 m ), relating the ob- 
served galaxy overdensity (6 ea i = «03i6 /"field - 1) to the mass 
overdensity and p = 3.5 x 10 10 Mq Mpc -3 for the cosmological 
parameters used in this paper. 

The weighted mean of the three density estimates in Sect. 
ED is H()3i6/"fieid = 3.3+° 4, giving an overdensity of 6 ga \ = 



2.3 + ° 5 r Taking V = 9.3 x 10 3 Mpc 3 (Sect.lO and b = 3 - 6 
dSteidel et all Il998l IShimasaku et al. , 120031) results in a mass 
for the protocluster within the observed volume of 4 - 6 x 10 14 
Mq. Because the size of the protocluster is unconstrained (e.g. 
Fig.1141. this mass estimate is a lower limit. 

However, the redshift range of protocluster galaxies is 
likely to be smaller than the redshifts for which the narrow- 
band filter is sensitive (see Fig. 1131 . Assuming that the three 
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outlying galaxies on the red side of the filter as field galaxies, 
the redshift range of the protocluster galaxies is 0.029 and the 
volume occupied by these emitters is 5.4 x 10 3 Mpc 3 . This es- 
timate of the volume does not take into accou nt the red shift 
space distortions caused by peculiar velocities JSteidel et all 
Il998l see below). Assuming that in total ~ 10% of the (can- 
didate) emitters are field galaxies (see Sect. 17.2k the density 
of emitters within this volume with respect to the field density 



is 1 + «5 gal = 3.3^5 x 0.9 x ||i 



5.1*? - ?. In this approach, 



the relation between the mass overdensity 5 m and the observed 
galaxy overdensity 5 ga \ is JSteidel et all fl998): 



l+W B = C(l+&»i), 



(22) 



where C takes int o account the redshift space distortions 
JSteidel et allll998h . Assuming that the structure is just break- 
ing away from the Hubble flow, C can be approximated by 

C= 1+/-/(1+<U 1/3 (23) 

JSteidel et allll998l) . with / the rate of growth of perturbations 
at the redshift of the protocluster JLahav et all Il99ll) . / not 
only depends on z, but also on Qm and Qa. In the cosmology 
adopted in this paper ( Qm = 0-3 and Q\ = 0.7), / is close to 1 
at high redshift (z > 2. lLahav et ailll99ll) . 

Again taking b = 3 - 6, the computed mass of the proto- 
cluster is > 3 - 5 x 10 14 Mq. As both the redshift range of 
the galaxies and the size of the protocluster could be larger 
than is observed (see Sects. IF!T1 and l8.2> . this estimate is again 
a lower limit. The computed mass roughly corresponds to the 
virial mass of the Virgo cluster (e.g. Fouau e^et alll200ll) . 

9. Nature of the Lyn emitters 

What can we deduce about the nature of the Lya emitters 
in our field? Four of the 16 emitters (25%) detected in the 
ACS image are unresolved and may be narrow -line (FWHM 
< 1000-2000 kms -1 , e.g. iBennert etafl l2002l) QSOs. The 
number density of faint (21 < m R < 25.5) QSOs near 
z = 3 in the LB G survey of Steidel and collaborators 
JSteidel et all l2003h is 250 Q SO deg~ 2 in a redshift interval 
Az - 0.6 (Steidel et al., 2002). They observe a ratio of narrow- 
line to broad-line QSOs of A^narrow)/7V(broad) = 1.2 ± 0.4 
JSteidel et all 120021) . Assuming an overdensity of galaxies in 
our field of no3i6/ n fieid = 4.0 (Sect. |7~TT >. the predicted num- 
ber of QSOs in field is ~ 1. Indeed, one (broad-line) QSO is 
found (emitter #2487). Extrapolating the faint e nd of the QSO 
luminosity function given by Hu nt et all l| 2004) to m R > 25.5, 
the number of QSOs near MRC 0316-257 with 25.5 < m R < 
28.5 is calculated to be < 1, making the identification of the 
Lya emitters with QSO unlikely. Another reason against the 
classification of the Lya emitters as QSOs is the continuum 
slop e of narrow-lin e QSOs in the LBG survey of /3 — -0.4 
(Steidel et al., 2002). This is much redder than the median of 
the unresolved emitters (J3 ~ -1.43). We therefore conclude 
that all the Lye emitters are star forming galaxies 3 

3 We can not entirely exclude that some of the emitters harbour a 
(concealed) AGN. A 40 ks exposure by the Chandra X-ray observa- 
tory is scheduled for early 2005, which should help to estimate the 
fraction of the Lya emitters that contains an AGN. 



The next question is how the properties of the Lya- emitters 
compare to those of the LBG population as a whole. 

• Continuum luminosity: Besides that of the radio galaxy 
and the broad-line QSO, the continuum of the emitters is 
faint, the brightest emitter having m R — 24.24 (using m R = 
-2.5 10 log(/ v (A rest = 1700 A)) - 48.59) an d the faintest 
m R > 28.45. Using m, = 24.53 at z = 3.13 dSteidel et all 
1999), this corresponds to luminosities ranging from 1.3 L* 
to < 0.03 L*. Roughly 90% of the emitters are fainter than 
m R =25.5, the spectr oscopic limit for LB Gs. A similar per- 
centage was found bv lFvnbo et all ( 120031) . 

• Size: A comparison of the sizes of Lya emitters with those 
of LBGs at z ~ 3 suggests that Lya emitters are generally 
smaller (with r/, < 1.5 kpc) than LBGs which have an aver- 
age size of r/, ~ 2.3 kpc. 

• Color: The colors of the confirmed emitters are, on average, 
very blue. The median UV continuum slope is f3 — -1.76, 
bluer than t he average slope of L BGs with Lya emission 
ifi ~ -1.09; Shapl ev et all l2003fc . A large fraction of the 
confirmed emitters (~ 2/3) have colors consistent with neg- 
ligible absorption and could be dust-free starburst galaxies. 



Summarizing these properties, the Lya emitters are on average 
fainter, bluer and smaller than z ~ 3 LBGs. 

Various models have been proposed to explain the prop- 
erties of Lya emit ters. Based on rest-frame optical photome- 
try of LBGs, IShaplev etaD J200ll) concluded that LBGs with 
Lya in emission are "old" (age > few xlO 8 yr), while "young" 
(age < 10 8 yr) LBGs have Lya in absorption. This could be 
explained by the young galaxies being dusty which caused 
the absorption of Lya photons, while the older galaxies are 
more quiescent with less dust and superwinds which allow 
Lya photons to escape. Other groups have suggested that 
strong Lya emitters are young star forming galaxies which 
was derived from the blue colors and high equivalent widths 
of the Lya emitt e rs (e.g. LF96. iRhoads & Malhotral 12001 ; 
iKeel et all l2002t iMalhotra & Rhoadst 120021 iTapkenetal 



2004). Trying to connect these o bservations, various au thors 
(e.g. iFriaca & Terlevichl 1 19991 iThommes & MeisenheimeA 
2004) discuss a model in which a forming galaxy has two Lya 
bright phases: an initial phase when the galaxy has started the 
very first period of star formation and is still nearly dust free. 
Due to supernova explosions, the interstellar medium will be 
enriched with metals, and dust will form in the galaxy. This 
dust will absorb the Lya photons, extinguishing the Lya emis- 
sion line. A second Lya bright phase occurs at a later time when 
strong galactic winds facilitate the escape of Lya emission. 

The observations of Lya emitters in our field support this 
second picture, and the Lya emitting galaxies might be young 
star forming galaxies in their first starburst phase. This should 
be confirmed by deep infrared observations. Modelling the 
spectral energy distribution of the Lya emitters from the UV to 
the rest-frame optical should allow the discrimination of Lya 
emitters being either young dust-free galaxies or more evolved 
star forming galaxies with an underlying old stellar population. 
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10. Implications of a protocluster at z = 3.13 

10.1. Star formation rate density 

The total UV star formation rate density (SFRD) of the con- 
firmed emitters (excluding the radio galaxy) within the volume 
of the narrow-band filter is < 0.0127 + 0.0003 Mq yr" 1 Mpc~ 3 . 
The SFRD derived from observations of LBGs at z ~ 3 with 
m K 5, 27 is 0.0184 + 0.0034 M yr-' Mpc 3 JSteidel et all 
Il999l) . Using the same magnitude limit, we find a SFRD of 
0.0109 + 0.0002 Mq yr 1 Mpc 3 . This is a lower limit, because 
it does not include a contribution from the radio galaxy, the 
emitters with only a limit on their UV SFR are ignored and no 
correction has been made for incompleteness. Assuming that 
only 20-25% of the star forming, UV bright gal axies at z ~ 3 
have a Lya line satisfying our selection criteria dSteidel et all 
l2000UShaplev et alll2003b . then the SFRD around MRC 03 16- 
257 is roughly > 2.4-3.0 times higher than in the field, in agree- 
ment with the number density of the Lya emitters. 

It should be noted that the total SFRD in the protocluster 
mig ht be much higher. Fro m rest-frame UV and optical col- 
ors, ISteidel et all dl999h and IShaolev et all feOOll) found that 
the UV continuum is on average attenuated by a factor of ~ 5. 
Also, very dusty, obscured gal axies could be missing. For ex- 
ample, iDe Breuck et all (2004) found an overdensity of bright 
sources at 1 .2 mm wavelength in the field of the protocluster 
near the radio galaxy TN J 1338-1942 at z = 4.1. These ob- 
jects could contribute substantially to the total SFR within the 
protocluster. 

10.2. Enrichment of the intracluster medium 

It is interesting to estimate the metal production in the proto- 
cluster surrounding MRC 0316-257. At redshift z ~ 0.3, the 
intracluster medium (ICM) i n clusters has a me tallicity of 0.2- 
0.3 Zq fe.g.lMushotzkv "& LoewensteinHl997l). showing little 
evolution up to z ~ 1.2 dTozzi et all l2Q03t [Hashimoto et all 
EoollMaughan et alll2004HRosati et alll2004>~ 

The extinction corrected total star formation rate den- 
sity of UV-brig ht star forming galaxies at z ~ 3.1 i s ~0. 1 3 
Mq yr 1 Mpc" 3 JSteidel et all Il 999t ICiavalisco et all bOOA . 
The protocluster has a volume of 5.4 x 10 3 Mpc 3 and a galaxy 
overdensity of «03i6/«fieid = 5.1 (Sect. IS31 . This gives a to- 
tal SFR in the protocluster of ~3580 Mq yr _1 , ignoring any 
contribution to the SFR from very dusty, obsc ured starform - 
ing galaxies. Taking an average yield of 0.02 llLia et all l2002). 
this means that ~72 Mq of metals are produced every year. 
Taking 4 x 10 14 Mq as the mass of the protocluster (Sect. 
1 8. 3t and assuming a baryon fraction of Q/,/Qm = 0.17 (e.g. 
Sper gel et all [2003 ) and assuming that the star formation rate 
in the protocluster is constant with time, then enough metals 
can be produced to enrich the baryons in the protocluster to 
0.2 Zq at z ~ 1. However, a large fraction (> 90%) of the 
produced metals must escape the galaxies in which they are 
formed. A possible mechanism to inje ct the metals into the 
ICM are supernova- driven outflows (e.g. Heck man et allfl995l 
lAguirre e t all l200 l | >. which are frequently seen i n z ~ 3 LBGs 
(e.g. Pett ini et all 12001: Adelbergeret al., 2003). This simple 



calculation shows that the star formation rate in the protoclus- 
ter is high enough to enrich the ICM to the observed value at 
lower redshifts. 

10.3. High redshift protoclusters associated with radio 
galaxies 

Based on the high volume density of Lya emitters near MRC 
0316-257, which is a factor of 3.3*2';, higher as compared to 
blank fields, and the small velocity distribution of the con- 
firmed emitters (FWHM of 1 5 1 km s~ 1 ) compared to the width 
of the narrow -band filter (FWHM ~ 3500 km s _1 ), we conclude 
that the Lya emitters are located in a protocluster of galaxies 
with an estimated mass of > 3-6 xlO 14 Mq. A likely scenario 
is that this protocluster will evolve and form a massive cluster 
of galaxies. The radio galaxy at the centre of the protocluster 
has the properties expected of the progenitor of a massive cD 
elliptical. The clumpy appearance of the radio galaxy in the 
ACS image could be explained as a merger of smaller subunits, 
and is very si milar to HST observa tions of other z ~ 2 - 3 radio 
galaxies (e.g.lPentericci et allfl999l). 

Based on the observations obtained in our VLT large pro- 
gram, the protocluster around MRC 0316-257 is not unique. 
We have found galaxy overdensities around all five radio galax- 
ies (with redshifts between z = 2.16 and z = 4.1) that 
were studied to a simi l ar depth as the 0316 fie ld ( Kurk e^ll 
200 0l IPentericci et all kOOOal IVenemans etail 120021 120031 
Kur k et all 120041) . Each of these radio galaxy field has at least 
20 confirmed protocluster members, the velocity dispersions 
of the protoclusters range from 300 km s" 1 to 1000 km s 1 and 
the associated masses are > 10 14 Mq (Venemans et al .,2003). 
At an even higher redshift, z = 5.2, we found a similar over- 
density of Lya emitters near the radio galaxy TN J0924-2201 
dVenemans et all 120041 . In a future paper, we will describe 
and compare the properties of these radio galaxy protoclusters 
(Venemans et al. in prep). 
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Fig. 8. VLT narrow-band and ACS images of Lya emitters 
with a clumpy counterpart in the ACS image. The cutouts of 
the narrow-band image are ~ 8" on the side. The ACS im- 
ages are zoomed in on the centre of the narrow-band image. 
The grayscale ranges from 0.5 to 5 times the rms background 
noise. The emission to the left of the radio galaxy in the VLT 
image (and in the top-left corner of the ACS image) are from a 
foreground galaxy at z ~ 0.87 (see also Sect.HJ. 



Fig. 9. VLT narrow -band and ACS images of Lya emitters 
that remained unresolved in the ACS image. The cutouts of the 
narrow-band image are ~ 8" on the side. The grayscale ranges 
from 0.5 to 5 times the rms background noise. 
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Fig. 11. The ratio of recovered size over input size as a function 
of magnitude in the ACS image. The fainter and/or larger the 
objects, the more the size is underestimated. 




10 



SFR Ly „ (M sun yr') 



RG 



V 



QSO 





-1000 



1000 

Av (km/s) 



2000 



3000 



Fig. 13. Velocity distribution of the confirmed Lya emitters. 
The bin size is 150 kms -1 . The median redshift of the con- 
firmed emitters z = 3.1313 is taken as the zero-point. The 
velocities of the radio galaxy and the QSO are indicated with 
arrows. The solid line represents the selection function of the 
narrow-band filter, normalized to the total number of emitters 
found. The hashed histogram represents the velocities of emit- 
ters with absorption in their line profile. The three objects in the 
red part of the filter (with a velocity greater than 1500 km s _1 ) 
might not be part of the protocluster. 



Fig. 12. SFR calculated from the UV continuum plotted against 
the SFR computed using the Lya flux for the confirmed Lya 
emitters. 
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Fig. 14. Spatial distribution of Lya emitters around the radio 
galaxy MRC 0316-257 (denoted by a square). The confirmed 
emitters in the protocluster are represented by circles (emitters 
with a redshift smaller than the median, z < 3.1313) and dia- 
monds (emitters with z > 3.1313). The triangles show the posi- 
tion of the three Lya emitting galaxies with a velocity > 1500 
kms -1 from the median velocity of the emitters. The sizes of 
the symbols are scaled according to the velocity offset from 
the median, with larger symbols representing emitters with a 
redshift closer to the median redshift. The pluses are objects 
satisfying the selection criteria (see Sect. 13. 3> . but are not (yet) 
confirmed. The quadrangle represents the outline of the ACS 
image. 
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Fig. 10. VLT narrow-band and ACS images of Lya emitters, which are resolved by the HST. The cutouts of the narrow-band 
image are ~ 8" on the side. The grayscale ranges from 0.5 to 5 times the rms background noise. 
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Fig. 15. Spectra of emitter #344 (left) and #995 (right). 
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Fig. 16. Spectra of emitter #1 147 (left) and #1203 (right). 
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Fig. 17. Spectra of emitter #1518 (left) and #1612 (right). 
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Fig. 18. Spectra of emitter #1710 (left) and #1867 (right). 
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Fig. 19. Spectra of emitter #2487 (left) and #3101 (right). 
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Fig. 20. Spectra of emitter #3388 (left) and HzRG (right). 



